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Spontaneously created vortex-antivortex pairs are the predominant source of flux noise in high- 
temperature superconductors. In principle, flux noise measurements allow to check theoretical 
predictions for both the distribution of vortex-pair sizes and for the vortex diffusivity. In this paper 
the flux-noise power spectrum is calculated for the highly anisotropic high-temperature supercon- 
ductor Bi2Sr2CaCu208+i, both for bulk crystals and for ultra-thin films. The spectrum is basically 
given by the Fourier transform of the temporal magnetic-field correlation function. We start from 
a Berezinskii-Kosterlitz-Thouless type theory and incorporate vortex diffusion, intra-pair vortex in- 
teraction, and annihilation of pairs by means of a Fokker-Planck equation to determine the noise 
CI spectrum below and above the superconducting transition temperature. We find white noise at low 

frequencies ui and a spectrum proportional to at high frequencies. The cross-over frequency 

between these regimes strongly depends on temperature. The results are compared with earlier 
Q ■ results of computer simulations. 
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I. INTRODUCTION 



It is well known that the layered structure of cuprate high-temperature superconductors (HTSC's) leads to en- 
' hanced two-dimensional fluctuations. These fluctuations are partly eye to spontaneously created pancake vortex 
pairs in the superconducting Cu02 layers. There are several attemptsEI to describe these vortices starting from the 
Berezinskii-Kosterlitz-Thouless (BKT) renormalization group theoryE These approaches differ in their predictions so 
O , that experiments are needed to decide between them. 

Most experiments designed to test the predictions of BKT-type theories indirectly measure the temperature de- 
pendence of the renormalized interaction- This quantity can be obtained from the exponent a of the non-linear 
J> • current- voltage characteristics V oc /"Ela A second approach is to measure the linear resistivity above T c , which 
is related to the superconducting correlation length. It has been shown, however, that the derivation of the resis- 
tivity within the framework, of BKT theory is at best only valid in a narrow temperature range, which is probably 
i^J ' inaccessible experimentally!]! Thus, most of our experimental knowledge about vortex pair fluctuations is based on 

measurements of the temperature dependence of just one quantity. Alternative approaches would be very welcome. 
\^ • Apart from the renormalized interaction and the correlation length, BKT-type theories also predict the tempe- 
\ rature- and size-dependent fugacity of pairs, and, consequently, the distribution of pair sizes and the total pair 
density — at least below the transition temperature. A generalized approachEl yields quantitative results even above 
T c . It takes care of the correct counting of overlapping vortex-antivortex pairs and takes local- field effects in the 
screened interaction into account. In this way terms of higher order in the vortex fugacity y are introduced into the 
I ' Kosterlitz recursion relations,0 facilitating a description of the vortex system even above T c . 

Only few experiments sensitive to the pair density have been performed, most of them on magnetic-flux noise. In 
the absence of an external magnetic field, the flux at the surface of an HTSC sample is due to the vortices in the bulk. 
£j \ This flux is noisy since these vortices perform a diffusive motion, carrying their magnetic field with them. Only few of 
these experiments have been done on HTSC's, mostly on bulk YBa2Cu306+a (Y-123).Ll However, BKT-type theories 
and hence the approach presented here are probably not applicable to Y-123 since its anisotropy is too small. There 
are also noise measurements on Josephson junction arrays These arrays are discrete systems with a relatively large 
!h ' lattice constant— for which the continuum approach presented below is not suitable. 

Rogers et aln perform experiments on very thin films of Bi-2212 in the absence of an external magnetic field. 
(Apparently experiments on bulk single crystals of Bi-2212 have not been performed yet.) The authors find a flux- 
noise spectrum following a w -3 / 2 law for frequencies uj > lu c , where the characteristic frequency ui c strongly increases 
with temperature. 

In this paper we determine the effect of vortex-pair fluctuations in both bulk HTSC's and ultra-thin films (containing 
one Cu02 layer) on flux noise. To be specific, we consider the highly anisotropic compound Bi2Sr2CaCu20s+,5 (Bi- 
2212) in a vanishing external magnetic field under the assumption that the superconducting layers are coupled only 
weakly so that the dynamics of the vortices in one layer is independent of that in the other layers. We further 
assume a large density of similar pinning centers. We will find that the spectral density of flux noise is governed 
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by the temporal magnetic-field correlation function. Interestingly_|-the same correlation function also governs the 
contribution of vortp&-pair fluctuations to nuclear-spin relaxation,E9tj albeit at much higher frequencies. 

Ambegaokar et alS3 employ a Fokker-Planck equation to obtain the linear response of a superfluid film to substrate 
oscillations. In this equation the authors include the interaction between vortices within the same pair. Similarly, we 
also start from a Fokker-Planck equation. However, we solve this equation to obtain the full space- and time-dependent 
vortex correlations needed for the calculation of the magnetic-field correlation function. 

A similar vortex system is studied by Houlrik et allI3 They derive a relation between the flux-noise power spectrum 
and the dissipation due to the vortices described by a dielectric constant e. This—relation is valid in the limiting 
case of a large pick-up coil, i.e., for the flux through a_large area. Houlrik et a/.L3 perform computer simulations 
on a generalized two-dimensional, discrete XY modeltU to obtain e. The mentioned relation is employed to get 
the noise spectrum. It falls off as 1/ui 2 at very high frequencies ui and shows a l/ijli 2 dependence for smaller u>. 
The lu~ 3 / 2 power law is in agreement with Minnhagen's phenomenological approach.ll3 In the present paper results 
for a continuous two-dimensional Coulomb gas model are obtained by direct calculation as opposed to simulations. 
Furthermore, we consider the opposite limiting case of a small pick-up coil. 

The present paper is organized as follows: In Secjn] we define the flux noise power spectrum and express it in terms 
of a magne tic- field correlation function. In Sec. Ill we present a model which enables us to calculate this function 



and in Sec. IV we discuss our results. 



II. THE FLUX NOISE POWER SPECTRUM 



We have the following setup in mind: The small pick-up coil of a SQUID (superconducting quantum interference 
device) magnetometer is placed at the surface of a large HTSC single crystal or of an extended one-unit-cell-thick 
film. For epitactically grown samples, the most natural way to mount the input coil is on a {001} (ab) plane, sensitive 
to the field perpendicular to the layers. In the following we restrict ourselves to this case. The flux signal is measured 
in the absence of any external field or driving force. The spectrum is then obtained by Fourier transformation. 

The flux-noise power spectrum is given by the Wiener-Khinchin thcorcmJl3 
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S<f,(u;) = - dt 0(iWO) cos tot. (1) 
n Jo 

Here, 4>{t) is the flux through the effective area A g of the input coil. If the diameter of the effective area is smaller 
than the typical length scale of magnetic field changes, A a &, the field is approximately uniform over the area A e g, and 
we can write 

2 f°° ~ Fi 
SVH = -A 2 cS dt k zz (t) cos ut= J- A 2 eS k zz (uj). (2) 

K Jo V 7T 

We have thus reduced the problem to the determination of the Fourier transform of the magnetic-field correlation 
function 



k zz (t) = ho, z (r,t)h , z (r,0), (3) 
where h, ltZ {v, t) is the z component of the total magnetic field at the point r in layer n at the time t. 



III. MODEL 



A. General considerations 



In this section we present a model for the time-dependent local magnetic field in a layered superconductor in the 
absence of an external field. Results for a single layer are then obtained by means of a straightforward generalization. 
We assume that the Josephson coupling between the superconducting layers can be neglected as far as the dynamics 
of pancake vortices is concerned. Then the local field is due to spontaneously created pancake vortices in the layers. 
We assume that there are N vortices and N antivortices in each layer at any time, thus neglecting fluctuations of the 
vortex number. This is justified since we are only interested in the thermodynamic limit. 

We decompose the vortex system into the smallest possible vortex-antivortex pairs, using the enumeration algorithm 
given by HalpcrinJiZI i.e., we count the vortex and the antivortex with the smallest separation as a pair and then repeat 
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this step for the remaining vortices and antivortices. Let H„(r) be the magnetic field of a single vortex situated at 
the origin is, layer zero measured at the point r in the n-th layer. Here, we only need the z component of the field. It 
is given bytS 



H ( F ) = 005 expj- Vr2+nVi 1 (A) 

where 4>q is the superconducting flux quantum and s is the interlayer spacing. This expression holds for an infinite 
stack of superconducting layers. The field differs from this result outside the crystal, where it is not screened. However, 
if the pick-up coil is placed close to the surface the difference should be negligible. The two-dimensional symmetric 
Fourier transform of Eq. (W) is 



(5) 



For now we only utilize the fact that the field of an antivortex is just the negative of H n (r). The fields of the vortices 
and antivortices are superposed to obtain the total magnetic field h„ (r, t) , which depends on time only through the 
positions of the vortices and antivortices. 

We are interested in the correlation function k zz (t) as given by Eq. (||). The total magnetic field in layer zero is 
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h„(r,()^^(H„[r-r_„, + (()] 

n— 1 U— 1 

-H n [r -!•_„,„_(*)]), (6) 

where r„ tU+ (t) (r„^_(i)) is the position of the vortex (antivortex) of the is-th pair in layer n at the time t. We now 
assume that inter-pair correlations are negligible as compared with intra-pair correlations. We keep the correlations 
between the fields of the vortex and the antivortex of the same pair, however. This approximation is justified if the 
typical pair size is small as compared with the average distance between neighboring pairs. Under the same condition 
the (extended) BKT theoryO is applicable. If we further assume diffusive dynamics we can write down the following 
ansatz for the correlation function:E£l 

M*) = ^E / d 2 r' + d 2 r'_d 2 r + d 2 r_ 

n=l J 

x [H n<z (r - r' + ) H n>z {v - r + ) 

-iT„ iZ (r-r' + )ff„, z (r-r_)] 
xP(r' + ,r'_;r + ,r_;t)/(r+-r_). (7) 

Here, /(r) is the normalized distribution function of the pair separation vector r = r + — r_. We obtain this function 
numerically from the extended BKT theory of Refs. ^| and 11. For our calculations we use an approximate form of 



/(r), which incorporates the essential physics, cf. Subsec. Ill C| . The diffusive motion of the pairs is described by the 



time-evolution kernel or diffusion function P: P(r' + , r'_; r + ,r^;t) d 2 r' + d 2 r'_ is the probability of finding the vortex of 
a given pair in the area element d 2 r' + about r' + and the antivortex of the same pair in d 2 r'_ about r'_ at the time 
t provided that the vortex was at r + and the antivortex at r at the time zero. The indices n and v of the vortex 
positions have been omitted in Eq. (Q) since we are dealing with one representative pair. 

B. The diffusion function 

It may be instructive to turn briefly to the case of unbound pairs. In this case the vortices diffuse independently 
and the diffusion function separates into a product of free diffusion functions for the two partners of the pair, 

1 / \r\ -r. 12 



P(r'r'_;r+,r_;t) = — — exp 
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where D is the diffusion constant of a free vortex. By rewriting Eq. ([|) in terms of center-of-mass and relative 
coordinates, we can see that the center of mass diffuses freely with the diffusion constant Dcm = D/2, whereas the 
separation vector diffuses with D Te \ = 2D. 

Now we wish to take two important effects into account, namely the interaction between the two partners of a 
pair and the recombination of pairs. The latter effect is expected to destroy the correlation on the time scale of the 
recombination time. The center of mass of the pair should perform a free diffusion. The time-evolution kernel can 
then be written as 

= 2^«K- E ir ! ) p " ,<r ' rit) ' (9) 

The task at hand is to determine the time evolution of the separation vector, P rc i(r', r; t). ■— ■ 

To obtain P le \ we solve a Fokker-Planck equation containing the intra-pair interaction V .t3 The vortex-antivortex 
interaction is given bya 

For most pairs the dielectric constant e(r) is close to unity! so that we may replace q 2 /e by q 2 and write V(r) rj 
q 2 ln(r/r"o). The error thereby incurred turns out to be small as compared with errors due to, e.g., the uncertainty 
of the diffusion constant. Our approximation is best justified for small pairs, for which the interaction is screened 
only weakly. For temperatures significantly above T c many large pairs with strongly screened interaction exist and 
the approximation breaks down, while the extended BKT theory also becomes invalid. 

If the mobility and diffusivity are isptropic and constant in space and time, the diffusion (Fokker-Planck) equation 
in the presence of a potential V readst£l 

dt 



AV^ + ^(VV) • VP rcl + D Yel AP ieh (11) 



where the mobility /i rc i is related to the diffusion constant through the Einstein relation /x re i = Dvei/ksT. The initial 
condition is P re i(r,ro; 0) = S(r — tq). Inserting the logarithmic potential we find 

dP i r 
' '"' - 2^ rcl <z 2 <5(r)p. ol + fi^q 2 — • VP reJ + A-ciAP ro i. (12) 



dt 



The first term on the right-hand side contains a delta function. This term yields a positive contribution to the time 
derivative only at r = 0. Therefore, it causes a (5-function term to appear in P at r = 0. Such a contribution does not 
affect P for r 7^ 0. Since "pairs" with r = are recombined and do not contribute to the magnetic field, we may omit 
the first term in Eqj-dl2j). After introduction of polar coordinates, the diffusion equation can be solved by means of 
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a separation ansatzjiil 



P re i(r', r; t) = ( — I exp 



4vrP> rc i< V r J V 4D Icl t 



DC 



x e M«>'-¥>)j ( ] ( I.:-;) 



where 



,= kBT ~ q \ (14) 
' 2k B T K ' 

and I a is a modified Bessel function. The full time-evolution kernel is obtained by inserting the solution for P re \ into 
Eq. §. 

For t = the diffusion function P ro i is normalized to unity by construction. At later times more and more weight 
is expected to accumulate in the (5-term at r = while the overall norm remains constant. The weight outside of the 
central singularity is obtained by integration over two-dimensional space, 
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where 

7(0,2;) = / dte-H 11 - 1 (16) 
Jo 

is the incomplete gamma function.0 

The expression ([l5]) indeed approaches unity for t — ► 0, but decreases monotonically with time and goes to zero 
for t — > 00. In particular, it behaves as ||Prel|| oc t 1 for large t (note that 7 < —3/2). In Fig. |l| the weight ||P TO i|| is 
depicted as a function of time for various temperatures. The time is given in units of r 2 /A/J, Te iq so that the curves 
are invariant under change of fi rc \. The mobility fi IC \ is kept constant. 

As shown in Fig. Yu, there is a plateau in ||P re i|| for small times and a sharp drop in the vicinity of an annihilation 
time r a = r 2 /Afi rc \q. This is the typical time the separation vector needs to diffuse from its initial value r to zero. 
The curve ||P ro i||(<) is smeared out at higher temperatures. If the separation vector assumes the value r' = 0, the 
pair is trapped by the singularity. Then the magnetic fields cancel exactly and the pair has annihilated. For low 
temperatures the pairs tend to creep "downhill" into the potential well until they annihilate after a time of the order 
of r a . At higher temperatures the diffusive motion is generally faster so that the first pairs rccombinc earlier, but 
many pairs first start to grow and recombine later. 

Note that pairs are created at the same rate as they are destroyed. However, newly created pairs do not contribute 
to the correlation function since their positions are not correlated with the pairs still existing or already destroyed. 



C. Distribution of pair sizes 



Apart from the diffusion function, we also need to know the distribution function of the separation vector, /(r), to 
calculate the correlation function. Unfortunately the pair size distribution is known only numerically. 

In the direction parallel to the layers the magnetic field of a vortex changes on a length scale given by the penetra- 
tion depth A a h. Thus, the fields of a vortex and an antivortex with a separation much smaller than A a {, almost cancel 
each other. These small pairs do not contribute significantly to the correlation function. We utilize this observation 
by approximating the pair size distribution by an analytical expression which becomes exact for large pairs. The pair 
size distribution is intimately related to the pair ftig-acity y 2 of BKT theory, /(r) = y 2 (r)/r A . The modified Kosterlitz 
recursion relations of the extended BKT theoryBO predict that y 2 and the renormalized interaction described by 
the stiffness constant K approach a finite, temperature-dependend fixed point y 2 (oo), K(oo) for large length scales. 
Hence, we can solve the recursion relations close to the appropriate fixed point to obtain the leading behavior of the 
fugacity, and thus of the pair size distribution /(r), at large length scales. We find that /(r) cx l/r 2 ^ +4 with 

( = -2 + ttK(oo) + 2Tr 2 y 2 (oo). (17) 

From Ref. |^ we see that the exponent £ vanishes for T > T c and is positive and, to leading order, proportional to 
y/T c - T below T c . Details may be found in Ref. |l[ 

A reasonable approximation for the pair size distribution function is 

This function shows the correct behavior for large r and does not introduce irrelevant problems at small r. Since 
BKT theory neglects pairs of size r < tq, they are not counted in the total density N/F. The correct normalized 
distribution then reads 



v 2C + 6 1 



2nr 2 tf[l-2/(2C + 6)]-¥[l-4/(2C + 6)] 
1 - (r/r ) 2 



1 - (r/V ) 2 «+ 6 ' 
where ^(x) = V (x)/T(x) is the digamma function@ 



(19) 
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D. Correlation functions 



Now we have all ingredients to calculate the correlation function k zz . Equation (0) can be rewritten as 



k zz {t) = ^LjrJ d 2 R'd 2 r'd 2 Rd 2 



H rh z ( R' + j Hn,z (R + 2 



-ffn.z I R' + — ) H n<z 



P^(R'-R;OP r ci(r',r;0/(r), 



(20) 



where is the free diffusion function of the center of mass. To make this expression tractable numerically, we have 
to analytically evaluate as many integrals as possible. As noted above we need the temporal Fourier transform of the 
correlation function. With P Te \ from Eq. (|l^) we get, as shown in Rcf. [Tl| , 



k zz {ui) 



2N 8V2^ 



F D 



rel 



f 00 poo 00 

I d 2 kY,\HnA^)\ 2 J drr l -<f{r) £ J„ 

n,— l ^ rn — 1 . m odd 



x / dr'r' x +iJ m \K- \ Rel 



V 7 2 +r. 




D Iel < V7 2 +™ 2 



/fc 2 w 



(21) 



where r< = mm(r, r 1 ) and r> = max(r, r'). Taking into account the special form of the vortex field as given by 
Eq. (|J), summing |H„ iZ | 2 over the layers, and performing the integral over the polar angle of k we get 



k zz (uj) 



2N 8V2tt 



l2 „2 r°° 



F Ad ^KJo V + X- 



dk k 



exp (2sJk 2 + \-* -1 



^° m=l, m odd ^ ' 



00 

< I dr'r ll+1 J n 
Iq 



kr' 
~~2 



Rel 




(22) 



To describe a single layer we just have to replace the sum over n by one term, say for n = 1. This simply leads to the 
replacement of l/(exp[2s(fc 2 + \~ b 2 ) 1/2 ] - 1) by 1/ exp[2s(fc 2 + A^ 2 ) 1 / 2 ] in Eq. @. 

Equation (|2|) suggests that u c ~ Z? rc i/4A 2 b is a characteristic frequency of the correlation function since to only 
appears in the expression fc 2 /4 + iuj/D rc \ and the characteristic value of k is l/X a b because of the exponential. In 
fact X a b is the largest length scale in the problem so that D rc \/AX 2 ab is the smallest frequency where we expect the 
spectrum to show any feature. 

Of the parameters appearing in the rates the numerical value of the diffusion constant D re \ = 2D is least well 
known. Here, we briefly discuss vortex diffusion and its relation to piiping. In the absence of pinning the friction 
coefficient r\ of a vortex can be obtained from Bardeen-Stephen theory,E3 r\ = 0Q/27rc 2 ^ 2 pO». To take the anisotropy 
into account, one replaces £ by the coherence length within the layers, £ a fc. We thus havetil T] a b — ^o/^^^abPn = (1 1 
with the effective mass ratio e 2 — m/M < 1. The mobility /1 of a vortex is then /i — l/rj a bd, where d is the thickness 
of the superconducting layers. The diffusion constant is obtained using the Einstein relation, 



D = nk B T = 



^c 2 e ab {T) P nk B T 
4> 2 d 



(23) 



If one employs the Bardeen-Stephen formula the diffusion constant in Bi-2212 turns out to be of the order of 1 cm 2 /s, 
which seems rather large. However, since Bardeen-Stephen theory neither takes into account the discreteness of the 
quasi-particle spectrum in the vortex cores nor the apparent d-wave symmetry of the energy g&n, it may well give 
incorrect results for HTSC's. Measurements of D for HTSC's do not present a consistent picturer 2 ] Diffusivities from 
10~ 4 cm 2 /s to 10 2 cm 2 /s have been reported. . 
A large density of similar pinning centers leads to a thermally activated behavior of the diffusion constant ,E2l 



D = D Q exp 



k B T 



(24) 



where E p is the typical depinning energy. Matters are complicated by the observation that the depinning energy 
depends on temperature. Rogers et alB find the following empirical relation for Bi-2212 films: 



G 



(25) 

with Ep/ks ~ 1200 K. Other experiments also support a large value of the activation energy.@ These results only 
hold on time scales longer than the typical depinning time. For shorter times the description by means of a diffusion 
equation breaks down and has to be replaced by a model explicitly incorporating discrete hopping. 



IV. RESULTS AND DISCUSSION 

From Eq. (^J) and the correlation function given by Eq. ( |22j ) we immediately obtain the noise spectrum S$ = 
y/2j^Al s ~k zz {uj). 

For the numerical evaluation of S$ we employ Monte-Carlo integration. For each set of parameters and each value 
of the sum index m — 1,3,... we have performed 3 to 40 Monte-Carlo runs with 5000 sample points each. We use 
the distribution of the results of the individual runs to estimate the numerical error. We find that the summands fall 
off quickly for increasing m so that the term for m = 5 is smaller than the error of the m = 1 term. Hence, terms for 
m = 7, 9, . . . are neglected. 

We first consider bulk Bi-2212. For the Ginzburg-Landau coherence length and the magnetic penetration depth 
we set £,ab{T = 0) s» 21.5 A and X a b(T = 0) ~ 2000 A. To obtain the lengths at a temperature T we employ the 
Ginzburg-Landau formula A a (,(T)/A a h(0) = £ a fc(^)/£ab(0) = Ak /T c0 /(T c o — T), where T c o is the mean-field transition 
temperature. The density of vortices and the parameter £ are determined from the extended two-dimensional BKT 
theory of Ref. ^. This is consistent sinoe.-we have neglected interlayer vortex correlations throughout this paper. As 
noted above, the extended BKT theoryErllll should be applicable even in a temperature interval above T c . For higher 
temperatures, however, any description in terms of vortex pairs fails and a picture of free vortices is more appropriate. 
In this case we expect the spectrum to fall off as I/uj 2 *— .The parameter 7 is given by Eq. (fl4|). For the coupling 
constant q 2 we make the standard linear approximation! 2 ^ q 2 = ?o^s(3nO — T), where g 2 can be obtained from the 
known values of k B T c /q 2 {T c ) « 0.2053, T c « 84.7 K, and T c0 « 86.8 KE3 

Since the diffusivity is not well known, we first display in arbitrary units at six different temperatures in the 
vicinity of T c as a function of lj/D tc \ in Fig. ^. Exemplary error bars from Monte-Carlo integration are also shown. 
Displayed in this way the curves do not depend on D re i. The units of are chosen in such a way that = 1 for 
w/D le \ = 10 3 cm~ 2 . The absolute value of the noise power is thus not comparable for different temperatures. Because 
of this choice of units, the pair density, which is a simple factor in S^, does not enter the calculation. The calculation 
thus becomes independent of D lc i and N/F, which are the two most uncertain quantities. 

The spectra show a cross-over from white noise at low frequencies to a drop at higher frequencies. The drop is 
found-to approach the power law cc w~ 3 / 2 (the dashed line in the figure). The same behavior is found by Houlrik 
et alH3 in their simulations, except at very high frequencies. A lo~ 2 drop in that regime, as seen by Houlrik et at, 
is not found. However, we cannot investigate higher frequencies since the numerical errors start to increase rapidly. 
Since higher frequencies correspond to shorter probed length scales the vortices should eventually act as free particles, 
leading to a u>~ 2 power law. 

The two frequency regimes of white noise and a w~ 3 / 2 power law are separated by a characteristic value of lo/ D Te \. 
Just below T c , this characteristic value strongly decreases with increasing temperature, whereas the temperature 
dependence is weaker for T > T c . Since the only quantity in our calculations that shows a similar behavior is the 
exponent Q of the distribution function /, the main source of the temperature dependence of lo c has to be C- A more 
rapid drop of /(r), corresponding to smaller average pair size, leads to shorter recombination times and thus to higher 
characteristic frequencies. In this way measurements of S^,(uj) probe the distribution of pair sizes. For T > T c the 
characteristic frequency is of the order of 10 7 cm~ 2 x D IC \, corresponding to a characteristic length 3 x 10 -4 cm, which 
is indeed of the order of X a b(T c ). (Note that X a b diverges at T c0 > T c .) 

To be able to compare the flux noise at different temperatures, we have to take the temperature dependence of both 
the prefactor of and the diffusion constant D Te \ into account. As a result we show the absolute noise power for 
bulk Bi-2212 at constant frequency as a function of temperature in Fig. ^. For any temperature, the value of uj/ D IC \ 
is fixed by the requirement that ui/D IC \ — 10 7 cm~ 2 at T = T c , together with the known temperature dependence 
of D Te \, cf. Eqs. ( p4| ) and (p5|). The noise strongly increases with temperature, which we mainly attribute to the 
temperature-dependence of the vortex density 2N/F. The density increases approximately exponentially as more and 
more pairs are thermally excited.u Additionally, there is a kink at T c , which should be the result of the kink in the 
exponent £. Since flux noise is dominated by large pairs, the increasing exponent £ below T c leads to an additional 
reduction of the noise. The characteristic form of S,p(T) shown in Fig. |^ can serve as an indication of a BKT-type 
transition. 
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With the diffusion constant D determined fcpm Bardeen-Stephen theory,E3 the characteristic frequency lies outside 
the experimentally accessible frequency range.0 We have argued above, however, that Bardeen-Stephen theory may be 
unapplicable to HTSC's. Turning the argument around, one could determine D le \ from experimental spectra. Hence, 
experiments on bulk samples close to T c are called for. Voss and ClarkeEH argue that a spectrum with S$ oc 1/w 3 / 2 
is expected for u> > 2D/A e s due to diffusion of vortices out of the sampled area. Since we consider the case of a 
very small pick-up coil, the cross-over to 1/oj 3 / 2 is expected to take place at rather high frequencies. Thus the high 
cross-over frequencies may be the result of our assumption of a small coil. 

We now turn to ultra-thin films of Bi-2212. We use the same parameters as for bulk Bi-2212 with the exception, 
of the BKT temperature, which we choose as T c « 31 K to allow comparison with the experiments by Rogers et al.U 
Figure |^ shows the flux noise spectrum of a film at three different temperatures. The units are chosen as before. 
Again, the spectra show white noise at low frequencies and a drop at higher frequencies. The spectrum does not 
follow a w -3 / 2 power law in the frequency range considered here. Weak convergence of Monte Carlo data precludes 
calculations for higher frequencies. However, we have no reason to doubt that uj~ z / 2 behavior is eventually reached. 
The qualitative shape of the spectra agrees with Ref. ^|. 

The cross-over frequency is again found to decrease with temperature, consistent with our picture of larger and 
larger pairs with increasing .-temperature, which take longer to recombine. This result is in contradiction to the 
experiments of Rogers et aln and the simulations of Houlrik et al£3, who find an increasing cross-over frequency. 
These are results for the opposite limiting case of a large pick-up loop, however. The origin of the discrepancy is not 
yet clear. Note that the simulations suggest a vanishing cross-over frequency at T = T c , which does not seem to be 
consistent with experiment. 

To conclude, we have used a model which is based on the assumption of diffusing vortex-antivortex pairs and 
incorporates intra-pair interaction and pair annihilation to obtain flux-noise spectra for Bi-2212 single crystals and 
films. The spectra show white noise up to a strongly temperature-dependent cross-over frequency and 1/w 3 / 2 noise 
above. As a function of temperature, the noise shows a distinct kink at the BKT temperature T c . We have shown that 
flux noise measurements can yield information about the size distribution of vortex-antivortex pairs and on vortex 
dynamics, and can be used as an additional probe for a BKT transition. 
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FIG. 1. The weight ||P r ci|| of the diffusion function outside the central singularity as a function of time for temperatures 
k B T/q 2 = 0.05, 0.1, 0.15, 0.2, 0.25 (the steepest curve corresponds to UbT /q 2 = 0.05). Time is measured in units of r 2 /4/x re ig 2 . 
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FIG. 2. Double-logarithmic plot of the flux-noise power spectrum for a c axis oriented Bi-2212 bulk single crystal at 
T = 84.2,84.5, 84.6, 84.7,85.2, 85.7 K (from top to bottom). The critical temperature is T c = 84.7 K. The dashed line 
corresponds to the power law S4, oc cu^ 3 ' 1 ' 2 . Exemplary error bars are also shown. 
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FIG. 3. The flux noise power for a c axis oriented Bi-2212 crystal as a function of temperature at fixed frequency. The error 
bars from Monte-Carlo integration are also shown. 
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FIG. 4. Log-log plot of the flux-noise power spectrum for a ultra-thin Bi-2212 film at T = 28, 31, 34 K (from top to bottom). 
The critical temperature if T c = 31 K. The dashed line corresponds to S4, oc lu~ 3 ^ 2 . Exemplary error bars are shown. 
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